Nearly-logarithmic decay of correlations in glass-forming liquids 
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Nearly-logarithmic decay of correlations, which was observed for several supercooled liquids in optical-Kerr- 
effect experiments [G. Hinze et al. Phys. Rev. Lett. 84, 2437 (2000), H. Gang et al. Phys. Rev. Lett. 90, 197401 
(2003)], is explained within the mode-coupling theory for ideal glass transitions as manifestation of the p-peak 
phenomenon. A schematic model, which describes the dynamics by only two correlators, one referring to 
density fluctuations and the other to the reorientational fluctuations of the molecules, yields for strong rotation- 
translation coupling response functions in agreement with those measured for benzophenone and Salol for the 
time interval extending from 2 picoseconds to about 20 and 200 nanoseconds, respectively. 

PACS numbers: 64.70.Pf, 61.20.Lc, 33.55.Fi, 61.25.Era 



Optical-KeiT-effect (OKE) spectroscopy is a powerful tech- 
nique for the study of the dynamics of supercooled liquids 
The experiment provides a response function Xa(0 for 
times t exceeding a fraction of a picosecond. The function 
Xa(0 is proportional to the negative time derivative of a cor- 
relator ~ (A*(r)A)/(|Ap}. Here, () denotes canonical 
averaging, and the probing variable A is the anisotropic part 
of the dielectric function. The instrumentation was improved 
recently by application of heterodyne detection. As a result, 
it was possible to measure XA{t) for times up to 500 nanosec- 
onds, i.e., glassy dynamics was documented for the enormous 
time interval starting at the end of the transient and extending 
over more than five orders of magnitude It was shown that 
the evolution of the glassy dynamics of Salol upon decreasing 
the temperature T can be interpreted by the universal formulas 
derived within the mode-coupling theory of ideal glass tran- 
sitions (MCT) | lj 2J. The fit values for the various parame- 
ters have been found to be consistent with those obtained by 
other light-scattering techniques However, the univer- 

sal formulas do not describe the data for all times outside the 
transient regime. Rather, if T decreases, there opens a time 
interval larger than two orders of magnitude, which precedes 
the interval of validity of the universal formulas. Here, the 
response for Salol follows closely a 1/f-law, i.e., the coiTe- 
lator exhibits nearly a logarithmic decay: (|)a(0 °= 
101 ■ This intriguing feature was observed also for some other 
van-der-Waals liquids, like benzophenone (BZP). This holds 
with the reservation that the heretofore unknown relaxation 
process can be described more adequately by a power law, 
<1>a(0 ^ f °^ —t''', albeit with a rather small exponent b' Ql- 
Gaining an understanding of the indicated findings is a chal- 
lenge to all theories aiming to unlock a comprehensive de- 
scription of liquids. In this Letter, it will be shown that the 
measured complex relaxation scenario 10 |5| is a generic, 
though not universal, implication of the standard MCT. 

The MCT was proposed as a mathematical model for glassy 
dynamics 01, whose fascinating features are obtained as im- 
plication of bifurcation points of nonlinear equations of mo- 
tion derived for the density fluctuations. The basic bifurcation 
is a fold singularity describing a transition from ergodic to 
non-ergodic behavior if the temperature decreases through a 



critical value Tc- Using £ {Tc — T)/Tc as a small parame- 
ter, the long-time behavior of the coiTelators can be evaluated 
by asymptotic expansions 171. The leading-order formulas es- 
tablish the universal features of the dynamics. Comparisons 
of experimental results with these formulas [2J and tests by 
molecular-dynamics simulation 01 have shown that MCT 
is a serious candidate for an explanation of glassy dynamics. 
Nearly-logarithmic decay of correlations was predicted within 
MCT for states close to bifurcations of the cusp type as can 
be inferred from Ref. fioil and the papers cited there. It was 
pointed out 12^ (JJ that the OKE data might be fitted by the 
universal results for relaxation near cusp bifurcations. How- 
ever, one has to vary at least two control parameters in order 
to approach such singularity. A fit of the OKE data, which 
only depend on the single control parameter T, would require 
implausible fine tuning of the coupling coefficients entering 
the MCT equations. The explanation proposed below is based 
on the existence of the simple liquid-glass-transition singular- 
ity. Our results are valid for response functions of variables A 
that couple strongly to density fluctuations so that a so-called 
p-peak can occur in the susceptibility spectra 111, 121. 

Let us consider an MCT model describing schematically 
the density fluctuations by a single correlator (|)(f). It obeys 
the Zwanzig-Mori equation of motion 

a,2(j)(f ) + vdMt) + £2>(r) + £2^ /' m{t - t')d,,(i?{t')dt' = , 

Jo 

(1) 

where £2 and v are frequencies parameterizing the short-time 
asymptote, (j)(f) = 1-1/2 {Q.t f + 1 /6 vQ.V + (t^). The 
kernel m(f ), which represents the interactions of the density 
fluctuations, is modeled as 



m(f) = Vl(^(f) + V2(^2(f)■ 



(2) 



Here, vi ^ and ^ are the coupling coefficients. The 
state of the system is specified by a point V — (vi,V2) in the 
VI-V2 plane. For small V, the model describes liquid states 
where (^{t ^ °o) =0. For large V, one gets a nontrivial long- 
time limit / = (|)(f °o), < / < 1. Parameter / quantifies 
the arrest of density fluctuations in the non-ergodic glass state. 
There are lines of liquid-glass transition points = (vjjVj), 
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each characterized by some number A,, 1 /2 ^ A, < 1 . This 
number determines the critical exponent a, < a < 1 /2, and 
the von Schweidler exponent b, < b ^ 1, hy the relation 
r(l - a) Vr(l - 2fl) = A, = r(l + by-/T{l + 2b). We con- 
sider the line, which can be parameterized by Vj — (2A — 
1 ) /X^ , v| = 1 /A^ , 0.5 < A < 1 . Crossing the line, the long- 
time limit of (|)(f) jumps from zero to the critical value = 
1 — A. The distance of the state V from V is specified by the 
separation parameter a ~ [(vi — Vj) + (v2 — V2)/'^]/'^(l —f^) , 
which is negative for liquid states and positive for glass states. 
This model is the simplest one reproducing all values for the 
anomalous exponents, which can occur in the general theory 
lUill . The parameter appears as plateau value of the (|)(f)- 
versus-logf curves for V near or as relative strength of the 
a-peak of the loss spectrum Q]. 

The probing-variable correlator obeys the same gen- 
eral equation of motion and initial decay as (|)(f ), but Q., v and 
m(f ) have to be replaced by the corresponding quantities Q.a, 
Va and mA{t). The model for the kernel is 

mA{t)^VAW)^A{t), (3) 

where va > quantifies the coupling of the probing variable 
to the density fluctuations. If v^i/' > 1, the long-time limit of 
(|)yi(f) jumps from zero to the critical value = 1 — 1/ (va/'^) 
as V crosses the transition line. Originally, this model was 
motivated for the density fluctuations of a tagged particle il4ll . 
Also, the MCT equations for reorientational dynamics of a lin- 
ear molecule suggest an expression like Eq. Q |]JJ. Within 
the microscopic theory, the coupling coefficients depend on 
the equilibrium structure functions, which depend smoothly 
on T. Therefore, v'l, V2, and va are smooth functions of 
the temperature. The specified model was applied repeatedly 
for the description of experiments, as can be inferred from 
Ref. 116,1 and the papers cited there. 

The long-time decay of the correlators at the transition 
point is given by [^{t)-f'-']/h = [(j)A(f) - /^l/ZiA = (foA)" + 
0{t^^). The time to depends on the transient dynamics. It 
is determined by matching the numerical solution of ^(t) for 
large times to the asymptotic formula. The amplitudes are 
h = X, hA= A/(va/'^)- In the limit of small |o|, there appears 
a large time interval where \<^{t) — /' | is small. Here, one gets 
in leading-order [(^(f) - /']//! = [(j)A(f) " /aI/^^a = Cag{t/t^). 
This is a scaling law with a correlation scale Co = 
and a time scale fa = fo/|c?|'''^". The function g{t) is deter- 
mined solely by A. It exhibits the critical power law for small 
rescaled times f, g (f <C 1 ) = 1 /f ", and von Schweidler's power 
law for large rescaled times in the liquid, g(f 3> 1) — Bf* |7]. 
The scaling law implies one for the response functions 

lit)lh = XA{t)/hA = s^k{t/ta) , (4) 

where % = c^/ta- The master function k{t) = —dg/dt inter- 
polates between the critical power law, k{t <C 1) = , 
and, in the liquid state, the von Schweidler power law, k{t ^ 
I) = -/7B/f 



The VH-light-scattering spectra of Salol f3\ or orthoter- 
phenyl L17J show that is about 0.9. This means there is 
strong arrest of the reorientational motion at the ideal glass- 
transition points in these van-der-Waals systems. Integrat- 
ing XA{t) over time, one can determine <^A{t) and read off 
the plateau. We find = 0.90 ± 0.05 from the T = 260K- 
data for BZP |5] and = 0.93 ± 0.03 from those for Salol at 
T = 257K 1 2], corroborating the preceding conclusion. There- 
fore, we focus on XA(f) for large translation-rotation coupling 

VA- 




FIG. 1: OKE response Xa(0 measured for BZP for temperatures 
r/K = 251, 260, 290, 320 (full lines from bottom to top) |5]. The 
dashed straight line has slope x = 0.80. The dotted lines exhibit re- 
sults A' ■ XA(f). with XA(f) calculated for the MCT model defined by 
Eqs. |l]OJ with SI = Q.a= 0.2v = 0.2va. The factors and £1 are 
chosen (from bottom to top) as 3.5, 6, 9.5, 8, and 1.67, 1.43, 2.0, 2.0 
THz, respectively. The left inset exhibits the variation of the coupling 
constants with changes of the temperature, and the right inset the one 
of the separation parameter o. The dashed line marked sc shows a 
scaling-law result, Eq. J4}, calculated for A = 0.70. The line in the 
lower left inset exhibits transition points V^; the point O = —0.003 
for A = 0.70, is marked by a cross. 

^ Figure n reproduces the OKE response functions for BZP 
|5]. A dashed straight line of slope x = 0.80 shows that 
the data for T = 251 K exhibit power-law decay x{t) °^ 
^ for 2 ps^ f ^ 20 ps. There is a von Schweidler-law- 
like variation: (|)(f) — const. —t^',b' = 1 — x = 0.20. In 
the double-logarithmic representation, the scaling-law result, 
Eq. 0, appears as interpolation between a straight line of 
slope ( 1 + fl) > 1 , and one of slope ( 1 — /?) < 1 . Changing the 
scales SahA and is equivalent to a translation of the curve. 
The dashed line marked sc shows an example calculated for 
A = 0.70 (a = 0.33, b = 0.64) shifted to match the data for 
r = 25 1 K. It provides a proper description of the observa- 
tions for 60 ps < f < 6 ns. However, since 1 + a > 1 > x, the 
scaling law cannot account for the data for f < 60 ps. This 
conclusion |5] is not altered by choosing A or t^, So differ- 
ently. The dots are the theoretical results for Xa(0 scaled by 
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amplitude factors and scales £1, which vary somewhat with 
T. Notice that the calculated functions exhibit transient os- 
cillations for times below and up to 2 ps, as do the measured 
curves. This shows that the parameters Q., v, Q.a, Va ^"e cho- 
sen reasonably. The strong variation of Xa(0 with changes of 
T is due to the changes of the coupling coefficients as docu- 
mented in the left inset. The state V — (vi,v'2) shifts towards 
the transition line if T decreases. An extrapolation to a tran- 
sition point for A- = 0.70, which is near the cross, is consis- 
tent with the fit values for vi and V2- Parameter o exhibits 
an almost linear temperature dependence, and extrapolation 
to o = suggests Tc « 235 K. Uncertainty estimates for X and 
Tc cannot be given yet, since we did not study in detail the 
possibility for data fits by other choices for vi, V2, and va- 




logjgt[ps] 

FIG. 2: Analog results as the ones in Fig.Qfor data measured for Sa- 
lol for T/K = 247, 257, 270, 300, 340 |2]. The dashed straight line 
has slope x=1.15. The MCT results are calculated for 7.9 THz = Q.= 
2Q.A = IOVa, V = 0, and factors N = 25, 37.3, 27.8, 17.9, 19, respec- 
tively. The open circles in the insets are the fit values for the other 
measured temperatures |2]; the corresponding response functions are 
not shown to avoid overloading. 

Figure|2]shows the OKE response for Salol for the highest 
and lowest temperature measured and three temperatures in 
between |2|. The data are normalized to unity for t — 0.01 ps. 
A dashed straight line of slope x — 1.15 shows that the data 
for T = 247 K can be described by 1 /f '^ for the interval 

2 ps < f < 100 ps. Equivalently, the correlator shows a decay 
similar to the critical one, (])(?) oc 1 /t" with a' =x—l = 0.15. 
For the T = 257 K data, a similar behavior is found with 
X fvQ, i.e., (|)(f) exhibits logarithmic decay |2]. The MCT fits 
have been evaluated for temperature-independent parameters 
£1, Q.A, V, Va. The normalization constants exhibit some T 
dependence. The path V = (vi,v'2) extrapolates to a critical 
point for X = 0.73 (a = 0.31, 0.59). The o-versus-r re- 
sults suggest Tc « 245 K. The numbers for X and Tc are consis- 
tent with those obtained previously Cliii . For the identi- 
fication of the ideal transition point it would be helpful to have 
data available also for T below the estimated Tc. The varia- 
tion of Va closely follows an Arrhenius law: va °^ exp(7b/r). 



where the fit value 7b = 1000 K is not unreasonable for a van- 
der-Waals liquid. 
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FIG. 3: Response functions for the model defined by Eqs. QjS} cal- 
culated with fl = O.A = v/5=Va/5 and va = 20 for three states spec- 
ified byX = 0.70, (vi,v2) = (v^,v5)(l +£), e = -10"", n = 2, 3, 4, 
yielding o = 3e/10. The results for x(0 are shifted upwards by 
3 decades. The dashed lines are the scaling-law results, Eq. J4}, 
where to = 0.755. The dotted line shows the leading correction to 
the scaling-law for short times. The dotted line marked cc exhibits 
the Cole-Cole response function Xa^(0- The insets demonstrate the 
superposition principle for the long-time dynamics. 

Figure shows the two response functions, x(f) = 
—d(^{t)/dt and Xa(0 = -d(^A(t)/dt, for three Hquid states 
near the transition point with X = 0.70 and the correspond- 
ing scaling-law results. There is a second scaling law, the 
superposition principle for the a process, dealing with the 
dynamics for times within the von Schweidler-law regime 
and longer: x(0 = X(r/ra)Ao, XA(f ) = XA{t/t'a)/C t'^ = 
fo/|c?P, Y = 1/ (2a) + 1 /{2b). The insets demonstrate this re- 
sult. The solution at V including the leading correction reads 
[<^A(f) - /a] = /M(foA)"[l +KA{tQ/t)"] and a corresponding 
formula with index A dropped. The number /Ta = — 1 .5 is 
evaluated from the coupling constants vi, V2, va; it is negative 
and large since is high |7|. The result is shown as dotted 
line for %a- The corresponding correction for x(0> ^ = 0.020, 
is so small that it cannot be made visible in Fig. |3l This ex- 
emplifies that the range of validity of the universal formulas 
can be quite different for different functions. The scaling-law 
formulas explain quantitatively the evolution of the glassy dy- 
namics as exhibited by %(f ) and this for the complete regime 
outside the transient, Q.t > 2. The state « = 2 is marked by 



4 



a cross in the left inset of Fig. [2 The states of interest for 
the data interpretation have a much larger separation from the 
critical point. One finds here, as in previous work J3, fl^ . 
that the scaling-law description can explain all the qualitative 
features of %(f), but that a quantitative description does not 
work for \a\ > 0.05. The same conclusions are valid for Xa (t), 
albeit only for times with tQ. ^ 3000. The size of the cor- 
rection term Ka is so large, that there appears a three-decade 
time interval outside the transient where the universal formu- 
las cannot account for the MCT solution. Within the interval 
2 ^ logjQ(i2f) ^ 3.5, the leading-order-correction formula de- 
scribes the results. But for 0.5 ^ logjQ(£lf ) ^ 2, even this for- 
mula is insufficient. Within the interval 0.5 ^ logio(nf ) ^ 3.5, 
nearly-logarithmic decay is exhibited by (|)A(f)- 

The equation of motion for (|)a(0 is equivalent to XaIw) — 
-£2^/[co(co + /Va) - ii^ + i^Ato »m(co)], where XA(tt>) and 
mA(tt)) are the Fourier transforms of Xa(0 and mA{t), re- 
spectively, for frequency co. For low-frequency phenom- 
ena, |co(co + /Va)| < this simplifies to Xa(co) = 1/[1 - 
CO mA(co)]. For f/fo :2> 1 and f/fa *C 1, one gets up to next-to- 
leading order: mA{t) = ,fm+hm{tQ/tY[\ +/r,„(fo/f)''], where 

fm = VAffl, h,„ = VA{fhA+fAh), K,nh,r, = Va [/'^/Ja^A + 

/^/zA' + /i/ja]. There is a trend that the negative contribution 
of the first term in the bracket is canceled by the positive one 
due to the last term. Hence, because the correction term Ka 
for ^A{t) is large, the correction term K,n = 0.16 for mA{t) 
is small. Thus, the critical law is a good approximation for 
mA{t). This leads to the Cole-Cole formula for the P-process: 
- X87[l + (-'«/cop)1, Xt) = +/'/>'a], cop = 
[(1 + fin)/h,nT{\ -fl)]'/"/fo- The Umit co/cop < 1 reproduces 
the universal critical decay: x{t) oc 1/f '+". The limit Co/cOp » 
1 leads to a von Schweidler-law-like decay: oc 1/f It 
depends on the relative size of t^^, (Op and ' , which part of 
the P-peak spectrum dominates the nearly-logarithmic decay, 
X{t)°^ 1+a. The dotted line in Fig.|3lmarked 

cc exhibits x'^{t) for our model. The preceding reasoning for 
using the leading asymptotic expansion for mA{t) rather than 
for (|)a(0 is valid also for the microscopic version of MCT. 
The schematic model used here is merely the simplest exam- 
ple illustrating our derivation of the p-peak phenomenon. 

We have shown that the evolution of glassy dynamics as 
measured by OKE spectroscopy for two liquids can be de- 
scribed reasonably by the solutions of a standard schematic 
MCT model, and this for time intervals extending up to five 
orders of magnitude. The long-time part of the response func- 
tions can be explained qualitatively by the known scaling 
laws reflecting leading-order asymptotic solutions of the MCT 



equations. The leading-order correction formulas explain that 
for strong probing-variable-environment coupling there ap- 
pears a large time interval adjacent to the one for the transient 
motion where the scaling laws cannot explain the data. The 
nearly-logarithmic decay within this interval is explained as 
p-peak dynamics and is fitted perfectly by the solutions of the 
schematic model. Our theory implies the prediction of nearly- 
logarithmic decay of reorientational correlations for liquids of 
molecules with large elongation. Molecular-dynamics sim- 
ulations can test this and provide the structure factors for a 
microscopic calculation within MCT. Such work could sub- 
stantiate or falsify the preceding explanation. It remains to be 
explained why the p-peak phenomenon has not been detected 
in earlier studies. 
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